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We consider the effects of the Faddeev-Popov determinant in the Coulomb gauge on the confine- 
ment properties of the QCD vacuum. We show that the the determinant is needed to regularize 
the otherwise divergent functional integrals near the Gribov horizon but still enables large field 
configurations to generate IR enhanced running coupling. The physical gluon propagator is found 
to be strongly suppressed in the IR consistent with expectations from lattice gauge calculations. 
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I. INTRODUCTION 

Quantitative understanding of confinement and more 
generally of the dynamics of gluons at low energies re- 
mains as the major challenge in QCD. In the past few 
years lattice simulations and phenomenological studies 
have provided new insights into the nature of the low en- 
ergy behavior of the gluon propagator and role of gluons 
in forming the hadronic spectrum Q, Q- Since gluons 
can only participate in strong interactions, spectroscopy 
of hadrons with excited gluonic modes is of crucial impor- 
tance for investigations of confinement. It has recently 
been shown that hybrid mesons with excited quark and 
gluon modes should have properties similar to that of 
ordinary hadronic resonances and thus gluonic cxcita- 
tions may appear in the meson spectrum $ g & Q 0. 
Searches for exotic mesons have produced a few tantaliz- 
ing candidates H, 0, 0, 0, 0, and new experiments 
planed for JLab and GSI in light and charm meson spec- 
troscopy, respectively, are expected to produce a map 
of gluonic excitations. In this paper we address gluon 
propagation in the QCD vacuum. This investigation was 
prompted by recent lattice results indicating that in both 
covariant and Coulomb gauges the low momentum glu- 
ons do not propagate |lltllallttllflll8| . This is precisely 
what one expects for physical degrees of freedom ejj. two 
transverse gluons in a physical gauge 0, 1^3, H3, l22| . 
In the four-dimensional Euclidean formulation of covari- 
ant gauge QCD, however, the lack of IR enhancement 
in the gluon propagator contradicts the naive expecta- 
tion that the color confining force could be simply related 
to the gluon propagator. A popular, phenomenological 
approach to gluon (and quark) low energy dynamics is 
based on a truncation of the self-consistent set of Dyson- 
Schwinger equations. In many such approaches the gluon 
propagator plays a central role in providing the effective 
interaction between quarks, for example, it is used to gen- 
erate dynamical chiral symmetry breaking p3 |. A soft 
gluon propagator implies that confinement has to be de- 
scribed by other means. For example it has been argued 
that in covariant gauges the Kugo-Ojima confinement cri- 
terion for absence of colored non-singlets in the physical 
spectrum can be satisfied with a soft gluon propagator 



and an enhanced ghost propagator [1|, |2J, |25| |26|, |27| . We 
will show that this also seems to be the case in Coulomb 
gauge formulation. 

In a covariant formulation one sacrifices positivity con- 
straints and the Fock space representation and intro- 
duces additional (ghost) non-physical degrees of freedom. 
Alternatively, by relaxing the requirement of manifest 
Lorentz covariance it is possible to eliminate all non- 
physical components and study confinement and other 
low energy phenomena within the framework of quantum 
mechanical wave functions. Such an approach has obvi- 
ous, important implications for quark model based phe- 
nomenology. Furthermore, at finite density it allows for 
well established, diagrammatic, many-body techniques 
to be used. A many-body approach has proven to be 
successful in treating a variety of low energy phenom- 
ena in QCD. For example the random phase approxima- 
tion, which is typically (e.g. for electron gas) relevant at 
high-densities, in QCD is also applicable at low-densities 
and may result in a self-consistent realization of confine- 
ment [22. l2cl l29l | . Due to the long-range nature of the 
confining interaction, at low-densities quasiparticle ex- 
citations have infinite energy which eliminates colored 
states from the physical spectrum. Due to the presence 
of bare quark-antiquark pairs near the Fermi-Dirac sur- 
face, the quasiparticle vacuum breaks chiral symmetry 
and leads to a non-vanishing scalar quark density. The 
collective excitations of this quark-antiquark plasma cor- 
respond to the Goldstone bosons [30l l3lll32L l3l I3I l35| . 

The picture described above relies on the existence of 
a long range, confined quark-quark interaction. Such an 
interaction is expected to arise from the Coulomb oper- 
ator which in the Coulomb gauge Hamiltonian describes 
direct interactions between (color) charge densities. Un- 
like QED, where this interaction is simply determined 
by the distance between sources, in QCD it is a compli- 
cated function of the transverse gluon field and can not 
be thought of as a simple potential i.e. of the Cornell 
type [36(. The conjecture that the Coulomb operator is 
related to the confining interaction is based on the ob- 
servation that it is positive definite and vanishes at the 
Gribov horizon. The Grivbov horizon defines the bound- 
ary of the gluon field domain. A number of approxima- 
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tions have been developed to calculate the expectation 
value of the Coulomb operator and to verify this con- 
jecture US @ iHil In the process it has been 
realized that the Gribov region still contains physically 
equivalent field configurations. To what extent the neces- 
sary identification of the wave functional at these gauge- 
equivalent points modifies the expectation value of the 
Coulomb operator remains an open issue 39] . 

The standard approach to a many body system is to 
introduce a physicaly motivated ansatz for the wave func- 
tional and to define approximations for the evaluation of 
expectation values. We have followed this approach in 
Ref. (2^ where we generated the confining interaction 
but a specific assumption on the behavior of the gluon 
propagator at low energies had to be imposed. In par- 
ticular, the gluon dispersion relation which follows from 
minimizing the vacuum expectation value was solved self- 
consistently together with the confining interaction but 
a particular boundary condition was chosen for the so- 
lutions. In absence of the Faddeev-Popov determinant 
this was necessary in order to obtain a nontrivial solu- 
tion to the coupled integral equations. In this paper we 
will show how the Faddeev-Popov determinant can be 
included, how it fixes the low energy gluon dispersion re- 
lation and brings it in a qualitative agreement with lattice 
results. 

The paper is organized as follows: in Sec. II we give 
a brief description of QCD in the Coulomb gauge. As 
discussed above the main novel feature of this approach 
is the inclusion of the Faddeev-Popov determinant. The 
details of the many-body formulation are given in Sec. III. 
Renormalization is discussed in Sec. IV and is followed by 
numerical results in Sec.V. 



II. QCD IN THE COULOMB GAUGE 

The QCD Coulomb gauge Hamiltonian, defined by V • 
A a (x) is given by , 

H = H [A Q (x), n Q (x)] = H Q + H qg + H g3 + H g , + H c , 

. (1) 

with, n°(x) being the canonical momentum satisfying, 

[n?(x),^(y)] = -zS ah S^(V)S 3 ^~ y), (2) 

<5J?(V) = S** — V'V^/V 2 and in the Shrodinger repre- 
sentation given by 11° (x) = — i5 /SA a (x). The five terms 
represent the kinetic energy, the quark-transverse gluon 
coupling, the magnetic 3— and 4— gluon couplings and 
the instantaneous Coulomb energy, respectively. In this 
paper we focus on the gluon sector and thus will ignore 
quark degrees of freedom. The gluon kinetic term is given 
by, 



H = ~ J d 3 x [j- l n a (x)JIL a (x) + (V • A a (x)) 2 ] , 



with 



J = Det (1 - A) 



Tr log(l-A) 



(4) 



being the Faddeev-Popov (FP) determinant. The matrix 
A, in the FP operator (1 — A) -1 , is given by, 



A 



x,a;y ,b 



gfacbA c (y)V y , 



(5) 



and in Eq. (0} the trace is over the spatial: x, y and color: 
a,b,c indices. The FP determinant is the Jacobian of co- 
ordinate transformation from the canonical coordinates 
of the Weyl gauge, A fi ' a (x) = (0, V a (x)) with kinetic en- 
ergy given by, 1/2 J d 3 x(— i<5/i5V a (x)) 2 , to the Coulomb 
gauge fields, A a (x) defined through a gauge map, 



V a (x) -> A a (x 



= u((j))A a u((l>) + u-^JVii^. 

(6) 

The dependence of the Hamiltonian and wave functionals 
on the iV 2 — 1 Euler angles </>(x) can be eliminated using 
the Gauss's law constraint and results in the Coulomb 
energy term [40|. 

H C = \f d 3 xd 3 y ^-V a (x)^[A] x , a;y , b p b ( y ), (7) 

with p a (x) being the color charge density, in the absence 
of quarks is given by, 



p a (x) = / abc n b (x)A c (x), 
and the Coulomb kernel K given by, 



^[A] x , a;yib = g 2 



(1-A)- 



1 

V 2 



(8) 



(9) 



x,a;y ,b 



(3) 



More details of te derivation of the Coulomb gauge can 
be found in Ref. HHH- 

Functional integrals in the Coulomb gauge are per- 
formed over the measure, H 7L ^ a ,idA l ' a (-x.)J . The Faddeev- 
Popov determinant results from the nonlinear field trans- 
formation given in Eq. 10 and reflects the complicated 
topology of the field space domain. Furthermore it is 
well known that the gauge condition, V • A a = is not 
a complete gauge fixing and thus the mapping V — ► A, <j> 
is not unique. The unique solution on a gauge orbit can 
be defined as the absolute minimum of the functional, 
I[A,g] = J dx(A Q (x)) 2 minimized over a f39|. At the 
minimum of I [A], V • A a = and the FP operator is 
positive. The space of the absolute minima defines the 
Fundamental Modular Region, A as shown in Fig. 1. The 
boundary of A is a set of gauge fields which lead to degen- 
erate absolute minima. The fundamental region resides 
inside the so called Gribov region, Q, corresponding to 
all minima of I[A] and thus also satisfying the transver- 
sality condition. The boundary <9£1 defines a set of con- 
figurations for which the gauge mapping is singular and 
J\A\ = 0. In what follows we will primarily study the 
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A=0 
x 




FIG. 1: A schematic representation of the field domain in the 
Coulomb gauge. The point A is on the common boundary of 
the fundamental modular region (6 A) and the Gribov region 
(SQ) and corresponds to a coordinate singularity, J[A\ =0 



role of configurations near singular boundary points on 
dfl. Since there exist configurations for which <9A and 
dft overlap, the point A a (x) = lies in A and both A 
and Q are convex, some fluctuations around the null field 
reach the coordinate singularity, 5£l without leaving A. 
Furthermore it has been pointed out |4l| that it is the 
common boundary points which dominate functional in- 
tegrals. Thus even if field configurations outside of FMR 
are included these may not lead to substantial errors. 

In the limit J = 1, the kinetic term describes a set 
of coupled harmonic oscillators and its (unnormalized) 
ground state |0) is given by 



(A|w ) = exp 



d 3 k 

(2nj~ 3 



u; (fc)A a (k)A Q (-k) (10) 



with u>o{k) = k = |k| being the free gluon energy and, 



A Q (k) 



dxe- lkx A Q (x), 



(11) 



represents the normal modes. The FP determinant leads 
to a suppression of the ground state wave functional near 
the Gribov horizon. This can be illustrated using an 
analogy between the Weyl and the Coulomb gauge kinetic 
terms and a harmonic oscillator in cartesian and spherical 
coordinates, respectively. In A-dimensions, the S'-wave 
harmonic oscillator radial wave function satisfies, 



1 N 

It, 



i 2 2 
+ UJ Xi 



R(r) 



7— 7— 

Or Or 



2 2 

uu r 



R{r) 



(12) 

Here the jacobian is given by J — r N 1 ~ exp(— AHog r) 
and it vanishes at the boundary r — * of the domain 
of r. The ground state wave function, R(r) is finite at 
that boundary, R(r) — exp(—r 2 uj 2 /2) but the radial wave 
function defined by u{r) — J x l 2 R{r\ vanishes as r — ► 0. 
The Hamiltonian can be redefined to absorb the jacobian, 



H -> H = J- 1/2 HJ 1/2 



(p 2 r + Lu 2 r 2 )+V c , (13) 



were, p r = —id/dr and the additional potential is given 
by V c =_J- 2 [ Pr ,J} 2 /^-J^[Pr, \Pr,J]]/2. The Hamil- 
tonian H is hermitian with respect to a flat measure in 
the radial direcction, 



/ DxR(x)HR(x) 



drJR{r)HR(r) = / dru(r)Hu(r) 

(14) 

In terms of H and u(r) one effectively recovers the simple 
harmonic motion in one dimension (modulo Vc), except 
for a boundary condition, u(r) — > at a point r = 
corresponding to the singularity of the coordinate trans- 
formation. 

The behavior of the wave function near the Gribov 
horizon is strongly correlated with confinement. In the 
Coulomb gauge, at the Gribov horizon, J — the 
Coulomb kernel diverges and this can be interpreted as 
a manifestation of confinement [22] . However, since the 
radial wave functional vanishes at singular points, <5f2, of 
the coordinate transformation, functional integrals over 
color singlet states are expected to be finite. In Ref. 
we have explored this scenario, but we have not ac- 
counted for the boundary condition on the ground state 
wave functional. In effect we used a gaussian ansatz for 
u and not for R, i.e. the wave functional was finite at 
the Gribov horizon. To ensure that functional integrals 
over the FP operator converge we had to choose a par- 
ticular condition on the parameters of the ground state 
wave functional. Summarizing, the FD determinant is a 
crucial element of the QCD Coulomb gauge dynamics as 
it surpresses the integrands in functional integrals near 
the Gribov horizon. 

As long as gauge fields are within the Gribov region, O, 
the Coulomb potential is positive and it is possible to use 
a variational approach. The simplest variational ansatz 
for which diagrammatic expansion is possible is given by 
the functional generalization of the harmonic oscillator 
ground state which is equivalent to the quasiparticle ap- 
proximation. This approximation can be systematically 
improved through a cluster expansion and excited states 
can also be studied ,35| . In this context one often uses the 
formalism of second quantization which is natural when 
dealing with gaussian integrals over polynomials (Wick 
theorem). In our case, however, before one gets to this 
stage one has to deal with non-polynomial operators, e.g. 
J[A\ or -ff[A], and thus it is simpler to procceed with the 
Schrodinger representation. 



III. THE QUASIPARTICLE SPECTRUM 

In the non-interacting case, H = Hq the perturbative 
vacuum of Eq. I|10[) minimizes the energy density of the 
system, i.e. 



= 



i) 



duo{k) 



4 



d 1 f d 3 q 

~ duj(k) 4 J (2tt) 3 

Here V is the total number of gluo 
V = 8 aa 6% f d 3 x = (JVg - 1) x 2 x 
the quasiparticle spectrum we will i 

I 



2(w|w) J J 

I 



B Q (x) = V x A a (x) + gf abc A b {x) x A c (x)/2, and 
I 



E C = / DA«j dxdye-f^^^-V [p^)JK[A] mbP \y)] e S*<*P±Wl-K (lg) 



(15) 



q degrees of freedom, 
volume. To describe 
ise the same gaussian 



variational ansatz. The vev of the Hamiltonian becomes, 



E(u) = (u\Huj}/(uj\uj) = E k {uj) + E c (u), (16) 



where 



The ground state normalization is given by, 



wo; 



f DA a J[A]e-f dkuWAaWA ' < -- u) = (J[A}) 

We use the (• • •) to represent functional integrals over 
the ground state ansatz functional with the flat measure. 
After integrating by parts the kinetic and Coulomb kernel 
contributions can be written as ( [dk] = eflc/(27r) 3 ), 

E K = -1- [ M7e -/*^)A"WA"(-k) 

2{w\ui) J 

J{dk}uj 2 (k)A a (k)A a (~k) + J dx{B a (x))'- 



(20) 



and 



e c = - J^ / da* ( nf[dk^ e -/ dkw(fe)Aa(k)A " ( - k) 

2{w\u) J J 
/ , a (k 1 ,k 2 )^[A] k ^, bP fo (k 3 ,k 4 ), 

(21) 

respectively. The charge density is now given by, 

p a {k % , k rf ) = f abc tu(k 3 )A b (k t ) ■ A c (k 3 ), (22) 
and the Coulomb kernel by, 

K[A]^ a , b - J dxdye l ^+^K[A] mb e^ +k ^. 

(23) 

Even though details of the boundary of the functional 
integrals are not known, the partial integration is pre- 
sumably justified since the integrand vanishes as A — > oo 



and at the boundary of the Gribov region, J — ► 0. Com- 
pared to the harmonic oscillator example discussed ear- 
lier, the partial integration combines contributions from 
Vc and p 2 . in Eq. (|13fl to the vev of the Hamiltonian 
and express them as coordinate space integrals over the 
gaussian wave function. In our case the complication 
in evaluating functional integrals over DA is due to the 
nonlinear dependence of J\A\ and K [A] on the FP op- 
erator, (1 — A) -1 = (1 — A[A]) _1 . These integrals are 
performed by expanding such functionals in powers of A, 
performing Gaussian integrals of over polynomials in A 
and approximating them by products of two-point cor- 
relations. This is will be illustrated in particular cases 
below. To simplify the notation, the triplet of indices 
representing momentum, color and spin will be denoted 
by greek letters, e.g. a = (k, a, i) and a doublet con- 
taining a momentum and a color index by a — (k, a). 
The summation convention will be used with upper and 
lower indices differing by a replacement k — » — k, e.g. 
A a _ A ia (k), 



^A a A a = A a A a = J2J2[ [dk]A ia (k)6%(k)A ia {-k). 

ot a ij 

.(24) 

In this notation A from the FP operator can be written 
as. 



= gfacb I [dk](27r) 3 <5 3 (p k q) AC( ^' tq = \% $ A\ 

(25) 
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with 



A 7/3 — A 



(p,o) 
(k,c,fc),(q,&) 



(27r) d (p-k-q) J acb . 

p 

(26) 



S<J ai $Ja n J— o 



7 (F 



2~J 



(29) 



Evaluation of functional integrals is simplified by intro- 
ducing corresponding generating functional, In this paper we are primarily concerned with the instan- 

/ taneous part of the transverse gluon propagator, 
UpdA p F[A] A ai ■ ■ ■ A a ™ e~~ ^ 1 



and since, 



we obtain, 

(A ai --- A a2 F[A]) = 



A 



J 

2ui 



(27) 



H« p = (A'*ApJ[A])/(J[A]) = ^f. (30) 



(28) 

where the last equality follows from translational invari- 
ance and color neutrality of the vacuum. In the approxi- 
mation J — one has £l(k) = u>(k) and one obtains the 
propagator used in Ref. [22]. From Eq. (|29|l it follows 
that, 



J S_ y ' • 



,2^7 4c 7 I^J 



A J 

A-\ 

2lo 



)/(J[A]) = jr- 

lLO n 



^ + 2^-^--^ {J 

0J a OJP J=0 



A J 

Ah — 

2uj 



)/{J[A] 



, (31) 



with the first term representing propagator the absence given by, 
of the FP determinant (J = 1) and the second term 



J 



SJ a SJi 3 



(J 



A-\ 

2lo 



)/(J[A}) = -< 



2tu a 



(i-xy 



J[A])/(J[A\) 



^(1-A)- 1 
2ujp 



J[A])/(J[A}), 



(32) 



where [A a ]^ 5 = d\^ a s /dA a = A^ Q ff . This relation can 
be represented through an infinite set of coupled integral 
Dyson equations containing all dressed vertices. As ar- 
gued in Ref. [I2,|23], however, vertex corrections give a 
finite and small modification and will be ignored. The 
dominant contributions in both the IR and the UV re- 
gions of the loop momentum integrals over instantaneous 
propagators come from diagrams with a maximal number 
of soft Coulomb lines and a maximal number of primitive 
self energy loops, respectively. The primitive self energy 
is shown in Fig. 2 and is given by 



^ 2uj. 



7 



7 



7 



2uj. 



7 



_ \"1 Xa _ t0 

7C A P ~ /3 i Q' 



(33) 



and explicitly, 



1% = 



g 2 N c 



dk] 



(i-(q-k) 2 ) 

2a,(|k|)(q-k)2- 



(34) 



This self-energy is UV divergent and has to be renor- 
malized. We will discuss renormalization in detail in the 
following section. To proceed we need to introduce the 
expectation value of the FP operator, 

d a - p = ([(1 - \)- l T ? J[A])l{J[A}) = 5 & - p d & . (35) 

A few lowest order diagram contributions to this vev are 
shown in Fig. 3. At the two-loop order the first and sec- 
ond diagram in the second line dominate in the IR and 
UV, respectively and are retained. In higher orders the 
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X 1 



a a 



P 



FIG. 2: Primitive self-energy, i.e. the lowest order correction 
to the FP operator, or the Coulomb line, 7°. The solid line 
represents the bare Coulomb potential, (1/p 2 ). 



dominant contribution comes from a series of rainbow- 
ladder diagrams obtained by summing the class of dia- 
grams generated by these two lowest order loop diagrams 
and results in the following approximation to the Dyson 
series, 

d* $ = 5 & - p d & = S a + ([A(l - \)fpJ)/{J) 

= 5& fi + E 4f[^dT B = 5 a - [1 + I a d & ] , 



(36) 



d = 5 




where 



(37) 



and explicitly, with d & — d(q), I & — I(q), is given by 
9 



d(q) 



1-5/(9)' 



(38) 



Using the same approximation (of ignoring vertex correc- 
FIG. 3: Rainbow-ladder approximation to the Dyson equa- tions), the Dyson equation for the instantaneous propa- 
tion for the Faddeev-Popov operator gator becomes 



<S 6 



(J 



A H 

2u 



)/{J[A]) = 



^d^Ld 

2ui a 2lo/3 



E 



d\ n d 

2uj a p 



X"d^d 

2LO/3 



(39) 



and is shown in Fig. 4. Since neutrality of the vacuum gluon propagator, l/2f2. In terms of this propagator, the 



implies 

[X a dX <rp^ = X< a s d s X s p ^ = 26 a pF a , (40) 



we finally obtain, 

where F a — F(q) 7 is explicitly given by 

N c 



(41) 



F(q) 



[dk] (1 ^ k ^ 2) d(lkl)d(|q-k|). (42) 



We can now return to the calculation of the vacuum 
expectation of the full Hamiltonian. By minimizing with 
respect to u> this determines u>, and from Eq. I|41|) . the 



kinetic vacuum expectation value, Ek is given by 



+ 2 V aaP V% s QJa gjfj dJ ^ dJs 



d d d d y 



A-\ 

2lo 



(43) 



The second term originates from the square of the mag- 
netic field, B a = [Vx] a 1 A''+V' 7 a pA a A 13 , and it is shown 
in Fig. 5, 



E K /V 



[dq] 



9 2 N C 
32 



uj 2 (q) + q 2 
2VL{q) 




FIG. 6: Low order contributions to the 2PI gluon-gluon 4- 
point function. 



FIG. 4: The rainbow-ladder approximation to the Dyson 
equation for the transverse gluon propagator 









FIG. 7: The vacuum expectation value of the Coulomb op- 
erator 



FIG. 5: The vacuum expectation value of the kinetic energy. 
The 4-point function contribution comes from the B 2 term 



This magnetic contribution involves a transverse gluon 
4-point function, which as discussed earlier, is approxi- 
mated by the product of two 2-point functions, i.e. the 
gluon-gluon scattering amplitude shown in Figs. 6 is not 
computed. The Coulomb energy vev is shown in Fig. 8 
and is given by, 



E c 



Id d d 9 [ V 


2 dJ a dJp dJ 1 dJs 



jfJn 



A 



J 
2uj 



)P~j8 



(45) 



with the charge density given by, Eq. I|22l) . p = p 1 = 

p^ a pA a A f} , and 



Ec/V = lr / [dk][dq r (q-k)i' {l <k ' q ' J ' 



(fi(|k|)-F(|k|)-fi(|q|)+F(|q|)) 2 



fi(|k|)fi(|q| 



(46) 



Finally the gap equation follows from, = d(Ex 
E c )/duj a , 



N cg 2 /■ rJ , 1 ( 3 -(q-k) 5 



Cl(k) 

N c / , rjlUl , , A f , 2 ,K(c l -k)(Q(k)-n(q)-F(k)+F(q))(Q(q)+n(k)-F(k)+F(q)) 



.|[dk](l + (q-k) 2 ): 



4 J JV VM ; ' (q-k) 2 

(47) 



with and / satisfying 

K(q) = f(q)d 2 (q), (48) , l-(q-k) 



/( g ) = l + iVo/w ^ q ^ k)2 /(l q -k]). (49) 
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It is also instructive to analyze the single quasiparticle The calculation is straightforward although more tedious 
dispersion relation, due to presence of up to three contractions corresponding 

the vev of six field operators. The result is 

EJ1 = (JA a HAf})/(J). (50) 



E(q) 



2fl(q) 



n 2 + F 2 (q)+q 2 



N c g 2 



[dk] 



(3-(q-k) 2 ) 
fi(fc) 



^£ fuun^n ,^ 2 ^(q ~ k) (n(fc) - Q(q) - F(k) + F(q))(Q(q) + Q(k) - F(k) + F(q)) + 2Q 2 ( q y 
- 4 J LflKj^ + W -Kj j (q _ k)2 n 



(51) 



After combining with the gap equation one obtains 
£(,) = S!(«) 



n(fe)(q-k) 2 
(52 

which is identical to the expression found in Ref. |22 
modulo replacement to — > fi. 



IV. RENORMALIZATION 

So far we have been ignoring potential UV diver- 
gences. These divergences should be removed by renor- 
malizing appropriate operators and the coupling con- 
stant, g. It turns out that all four equations of interest, 
Eqs. P5|l. gTJ, |@7J|, (UjJ, require renormalization. These 
equations have to be regularized first and this can be done 
by cutting-off the momentum integrals, / [dk] — > / [dk]. 
The physical, renormalized solutions, d(q),f(g), oj(q) and 
Q(q) should be A independent. We will first discuss 
renormalization of the expectation value of the FP op- 
erator, d(q). Assuming that a renormalized solution for 
£l{q) has been found, the equation for d(k) is renormal- 
ized by adjusting the bare coupling, g — > g(A), i.e. the 
renormalized vev of the FP operator, d(k) will play the 
role of the running coupling. The A dependence of <?(A) 
is determined by the UV behavior of Eq. I|38[) , 



rfg(A) _ (3 g 2 (A)d(A) 
dA (4tt) 2 n(A) : 



(53) 



with (3 = SNc/i- In this and all other renormalization 
group equations we keep only relevant and marginal con- 
tributions, i.e. no power corrections, 0(p n /A n ) with 
n > are included since they do not require renor- 
malization. Since (3 > and physical solutions re- 
quire d(k),fl(k) > the solution of Eq. (|53|l vanishes 
in the limit A — > oo. In the limit k — > oo the integral 
I(k = A, A) given by Eq. I|38|l . with the second argument 
referring to the upper limit of integration, is finite. Thus 



in leading logarithmic approximation, d(A) — > g(A) as 
A — > oo. Furthermore, from Eqs. I|41() it follows that for 
large q, q ~ A, 



dtt(q) diu(q) 
dq dq 



0{d 2 {q)) 



duj(q) 
dq 



+ 0{g 2 ). (54) 



Similarly from Eq. (|47|l . to leading logarithmic approxi- 
mation, we find du>(q)/dq = 1 + 0(g 2 ). Thus finally, 



A 



dg(A) 



T g 3 (A) + 0(g 5 (A)), (55) 



dA (4tt) 2 - 
which, ignoring the terms 0(<? 5 ) has a solution given by 



3(A) 



s(m) 



(l+(4^S 2 Gu)log(A 2 /M 2 ) 



1/2' 



(56) 



The asymptotic behavior as A — > oo is therefore given 
by, 



5(A)- 



47T 



/3V2 log 1/2 (A 2 ) 



(57) 



The renormalized equation for d{k) is completely speci- 
fied once g(fi), the value of the coupling at an arbitrar- 
ily chosen renormalization scale, /i is fixed. It should 
be stressed, however, that this solution is valid only to 
within terms of the order of 1/ log 3 / 2 (A 2 //i 2 ). In practi- 
cal applications we will be renormalizing at a low energy 
scale, fi e.g. related to the string tension or the glue- 
ball mass and thus such corrections become unimportant 
as A — > oo. For relevant operators, however, as we will 
see below, such logarithmic corrections are multiplied by 
positive powers of A and thus cannot be neglected. 

For practical (numerical) applications we have found a 
different, momentum subtraction renormalization (MSR) 
scheme to be more practical. In this scheme the renor- 
malized equation for d(q) is obtained by subtracting 
Eq. at q = /Lt, 
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1 



1 







d(q) d(lJ,) (Any 



(k-q) 

l Jo 



^ tf 3_(l-(k.qf)d(q-k) 



Q(k) (q-k) 



+ (q -> a*)- 



(58) 



In this renormalization scheme the coupling constant is 
therefore given by 



1 



1 



9ms r(A) d (v) 



(4tt) 2 



tffcfc 



2 3(l-(k-/i) 2 )d0u-k) 



(M-k)a' 
(59) 



and g(A) = g l jvfSi?(A) to within corrections of the order 
0(/x/A) i.e. they agree asymptotically. From now on we 
will drop the MSR subscript. Finally for fc ~ A — > oo 



d(fe)= 5 (A) [l + 0(g 2 (A)log(A 2 /fc 2 ))] 



(60) 



which is expected, as discussed above. 

We now proceed to discuss renormalization of the equa- 
tion for f(k). Physically f(k) represents additional con- 
tributions to the vev of the square of the FP operator, 
((1 — A)~ 2 ) ~ d 2 f not present in the expectation value, 
((1 — A))~ 2 ~ d 2 . Any UV divergent contribution to 
/ should therefore be renormalized by renormalizing the 
operator g 2 /(l — A) 2 , since the operator g/(l — A), has al- 
ready being renormalized. This is done by a renormaliza- 
tion constant, Z%(A) introduced by replacing the com- 
posite Coulomb operator, K[A) — > Zk(A)K[A]. Using 
the renormalized Coulomb operator the equation Eq. I)49|) 
for f(k) becomes, 



/(A;) = Z K (A) 



(4tt) 5 



J7>2 3 (l-(k.q) 2 )d 2 (q-k)/(q-k) 

(k - q) / dkk 4^m w^w — • 



(61) 



r 



and in the limit A —> oo one obtains, 



A 



dZ K (A) 
dA 



P rf 2 (A)/(A) 
(4tt) 2 0(A) 



Zk(A), (62) 



which in the leading logarithmic approximation has a so- 
lution given by 



Z K {A) 



log 1 / 2 (A 2 / / z 2 )' 



(63) 



Choosing a value for Z(/j,) at some UV point fixes the 
renormalized equation for f(k). As in the case of the FP 
determinant, we will employ the momentum subtraction 
renormalization scheme, which leads to, 



/(*) - m 



p r 1 



(4tt) 5 



A 



(k-q) / dkk 1 - 



3(l-(k-q) 2 )d 2 (q-k)/(q-k) 



4 n(k) 
i 



(q-k) 2 



(64) 



resulting, in the MSR scheme in Zk(A) given by, 



Z K (A) = f(p) - 



(4tt 



2 , (k. A )/' A ^?- (1 - (fc - A)2)d2(/i - k)/(/i - k) 



4 Q(k) 
I 



(M - k) 2 



(65) 



As expected, for UV values of k ~ A — > oo we find Much of the discussion on renormalization of d and / has 

f(k) = Z K (A) [l + O (g 2 (A) log(A 2 /fc 2 ))] . (66) 
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already being given in Ref. [22j . It should be noted that 
if for large k, d 2 (k)/£l(k) < l/log(fc)™ with n > 1 there is 
no renormalization for f(k). Our analysis suggests that 
n = 1 thus the unrenormalized equation for / has a sub- 
leading log(log(fc)). We suspect that this is an artifact of 
the rainbow-ladder truncation. 

As long as one works with the leading logarithmic ap- 
proximation, n(k) = u>(k) = fc, and there is no effect of 
the FP determinant on d or /. The inclusion of the FP 
determinant has an effect on the low momentum behav- 
ior of d and /, but it also introduces a new divergent 
integral, F{q) in Eq. I|42l) . Since the origin of F is the 
FP determinant J : it is the FP determinant that has 
to be renormalizcd in order to make Q(k) finite. The 
renormalized FP determinant should by chosen as, 



J 



Je 



E 



8u a A a A a + - 



(67) 



Here • • • stands for higher powers of the field operators, 
however, since including the FP determinant within a 
gaussian approximation to the functional integrals only 
the quadratic term needs to be retained. It can be easily 
verified that replacing J by Eq. (|fi7|) leads to the replace- 
ment, 



F(q)->F(q,A)+5oj(A), 



(68) 



where F(q, A) stands for the integral in Eq. I|42|l with 
the upper limit set to A. The counterterm Su>(A) will 
be chosen to make Cl(q) UV finite. Since F(q) has mass 
dimension of one, in general one expects two counter- 
terms would be needed, one proportional to A and the 
other to one power of the momentum. From the UV 
behavior of the integrand in Eq. I|42[) it follows, however, 
that only the first is needed and we obtain, 



dSm(A) 
dA 



whose solution is given by 



5u)(A) = 8uj{fj) 







(4tt) ; 



-rf 2 (A), 



(69) 







(4tt) 2 



dkd 2 {k). 



(70) 



We note here that corrections to the leading asymptotic 
behavior d{k) ~ g(k) cannot be neglected here since for 
F(q) they result in terms of 0(A). Thus it is necessary 
to keep d(k) rather then g in the renormalized expression 
for F. As in the case of d and / in the following we will 
use the MSR scheme for Sl(q) which gives, 



fl(q) — — uj(q) + u(n) 







(4tt) 2 



(k-q) 



dkk 



2 3(1 



(k-q) 2 ) 



4 (q - k) 2 



d(k)d(q-k)-(q^ v), 



(71) 



r 



with asymptotic behavior, k ~ A — * oo given by, 



fi(fc) = u(k) + g 2 (A)A [l + O (g 2 (A) log(A 2 /fc 2 ))] , 

(72) 



and with 5u>(A) in MSR given by, 

Slu(A) = - u(n) 



(47T) 



(k-q) 



dkk 2 - 



2 3(l-(k-/i) 2 ) 



4 (p-k) 2 



d(k)d(p - k). 

(73) 



The gap equation is the one which cannot be renormal- 
ized in a simple way like the previous equations. This is 
due to inconsistencies in the approximation used. Specif- 



ically, the gap equation is derived by taking the func- 
tional derivative of the energy expectation value with re- 
spect to a;. In the second integral in Eq. l|47Jl we have 
retained only the derivative of the gluon lines and not 
of the Coulomb kernel. The former leads to terms in 
the integrand proportional to F thus formally of 0(g A ). 
Similarly derivatives of the Coulomb operator d 2 f lead 
to terms proportional to d 4 //fi 2 i.e. also of 0(g 4 ). Thus 
if terms proportional to the difference f2 — u> are kept in 
the numerator of the gap equation it would be necessary 
to include derivatives of the Coulomb kernel. However 
all these 0(g A ) terms involve two-loop integrals and for 
simplicity will be neglected. The simplified gap equation 
then reads, 



\q) = q 2 + 2F(q)u(q) 



J 



N c g 2 



[dk] 



(3-(q-k) 2 ) , N C 



oj(k) 



-± /[dk](l + (q-k)< 



A-(q-k) u 2 {k)-L0 2 {q) 
(q-k) 2 Lu{k) 



(74) 



i.e. is the same as in Ref. [2^ except for the term in- volving F in the r.h.s and all terms are of 0{g 2 ). This 
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simplification is justifiable since our goal is to study the 
effect of the FP determinant on the low momentum prop- 
erties and thus possible modifications of UV behavior are 
largely irrelevant. 

In a covariant formulation the renormalized theory has 
the same operator structure as the bare one. This is not 
the case in the Hamiltonian approach. Renormalization 
introduces non-canonical operators. The strength of such 
operators is determined by the cutoff. The gap equation 
has a quadratic divergence which is to be renormalized 
by a gluon "mass" counter-term in the Hamiltonian, 



1 



SH(A) = -m 2 (A)A a A c 



(75) 



This is the only relevant operator e.g. of dimension two. 
The constant m 2 (A) is fixed by requiring that the gap 
equation leads to a A-independent solution. Thus we in- 
sist that uj(k) is A-independent and this guarantees that 



any divergence of an operator matrix element calculated 
with respect to the state \oj) will be associated with the 
operator itself and not the sate. The counter-term SH(A) 
contributes to the r.h.s of Eq. I|74|l with m 2 (A) and from 
the UV behavior of Eq. we find, 



dm 2 (A) 



(3 



dA (4tt) 2 
whose solution is 

m 2 (A) — m 2 (/i) 

/3 f A dk 



2 5 2 (A) 



A 2 



d 2 (A)/(A)u;(A) 



(76) 



(47T) 2 A w{k) 



[2g 2 (A)k 2 +d 2 (k)f(k)uj 2 (k)] 



(77) 



In the MSR scheme the gap equation then becomes, 



J 



\q)-L0 2 {p,)-q 2 -F{q)L0(q) + F^UJ^) + fi 2 = 



13 f 1 



(47T) 



(k-q) / dkk 2 [uo 2 {k) -u 2 {q)] 



3 (l + (k-q) 2 )g(q-k) 
8 w(k) (q-k) 2 



(78) 



with the asymptotic behavior for k ~ A — > oo given by, and the counter-term m 2 (A) in the MSR is given by, 



u 2 (k)=q 2 [l + 0(g 2 (A)\og(A 2 /k 2 ))] 
' 13 g 2 (A)Z K (A)A 2 , 



(4tt) 



(79) 



m 2 (A) = uj 2 (fi) + fi 2 



1 



-2-/^ 5 2 (A) / dkk' , 
(4tt) 2V J uj(k) (4tt) 2 



(k- A) / dkk 2 W{k)- uj 2 (p)} - 
l Jo 8 



3(i + (k-A) 2 )^-k) 



w(k) 



(/x-k) 2 



(80) 



The renormalized equations for the vev of the FP opera- 
tor, the corrections to d 2 needed to obtain the Coulomb 
potential, the gluon propagator and the ground state 
wave function are given by Eqs. (JjJBJ, dJ> an d GEl, 
respectively. These equations depend on four parame- 
ters, the renormalization constants, <i(/z), /(/i), Q{fJ-) and 
w(/i). In the following section we will study the solutions 
of these equations and their physical interpretation. 



V. RESULTS 

As discussed above there are four constants which need 
to be fixed. This can be done, for example, by comparing 
the Coulomb potential in position space, 



d 3 k 



tk . x rf 2 (k)/(k) 

k 2 



(81) 



with the lattice, static quark-antiquark potential. This 
procedure was used in Ref. p3 |. Unfortunately, the de- 
pendence of V e ff on the renormalization constants is 
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complicated thus fitting the lattice potential will not 
necessarily provide much physical insight. Furthermore 
it has recently being shown that there are differences 
between the lattice Coulomb potential and the Wilson 
loop potential [42|. We will thus proceed by simplifying 
the resulting equations and imposing constraints on the 
renormalization constants. The main difference between 
present analysis and what was done in Ref. [22j has to do 
with inclusion of the Faddeev-Popov determinant. Our 
goal here is to investigate the role of the FP determinant 
through the behavior of fi(fc) at low momenta. If the 
determinant is omitted, one has Q(k) = u)(k) and in this 
case the remaining three equations, for d(q), Eq. (J5SJ, 
f(q), Eq. and uj(q), Eq. (|7H|l . were analyzed in 

Ref. |22 • These equations have solutions provided uj(k) 
is finite as k — > 0. If ui(k) — > then the equation for 
d(k) will develop a pole at a finite, positive value of mo- 



J 



mentum and if u>(k) — > 00 as k — > then for a confining 
potential, K(k) — > l/k a , with a > 2 the gap equation 
has no solution. A renormalization condition at fi = 0, 



w(0) 



was therefore imposed with m g fixed by the 



Wilson loop string tension. A simplified set of equations 
can be obtained by making an angular approximation, 



|q - k| -> q6{q - k) + k,9{k - q) 



(82) 



where q = |q| and k = |k|. In Ref. [22j it was shown 
that the angular approximation lead to results which are 
very close to the exact numerical solutions. We will thus 
follow this approximation here since it allows us to con- 
siderably simply the numerical analysis. Using the angu- 
lar approximation the equation for the FP operator, d(k) 
becomes, 



1 1 



(3 f%,pd(q) , paw , a r^iw 



dk- 



dk- 



dk- 



d{q) (4ir)*J q *n(k) (4tt) 2 J„ n(A) (4tt)2 7 ^(fc)- 

I 



(83) 



the gluon propagator function, fl is given by 



n( q ) = no*) + w(flf) - + -/-^ /' dk%d(k)d( q ) + -JL; r dkd 2 (k) - JL; T dk^d(k)d^), (84) 

(4tt)^ 7 (4tt) 2 y g (4tt) 2 7 ^ 



and the gap equation for oj(q) becomes, 



uj 2 (q)-uj 2 (fi)-q 2 - F(q)ui(q) + F(p,)u>(fi) + fJ 2 = 



(4tt) 2 J q 2 Oj(k) 



dkK{k) 



L0 2 (k)-L0 2 (q) 

Uj{k) 



(85) 



r 



As discussed above if one ignores the FP determinant it is 
not possible to choose an arbitrary renormalization con- 
dition for w(0). Furthermore in this case there is a critical 
(maximum) value of d(/i) — d c (fi) = 47r^/3/3cij(0)/5/i for 
which Eq. (|83|l has a solution. Appearance of such a criti- 
cal coupling is an artifact of the approximation (rainbow- 
ladder truncation) used in evaluation of the expectation 
value of the Hamiltonian. The problem can be illustrated 
by considering the following integral, as a schematic rep- 



resentation of the function integrals representing the vev 
of the FP operator, 



1(g) = / dxj{x) : 



gx 



Here x represents the gauge potential, J(x) 



(86) 



Jog(l-x) 



plays the role of the FP determinant, and 1/(1 — gx) the 
FP operator. It one sets J = 1 the integral becomes 
divergent at x = 1 unless g = 0. The rainbow-latter 
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FIG. 8: Plot of the 

(r.h.s. of Eq. @2))/fi(0) 

31/25 - 25/125^ + 144(1 - 
z = d 2 (fi)fi//3n{0). 



function f(z 
-31n(l + 

1A/ 



5s/3)/5 + 

52/3) /625z 2 , with 



FIG. 9: A plot of the function f(z) = 1 
(r.h.s. of Eq. (ED)/n(0) = -31n(l + 5z/3)/5 + 31/25 
25/1252 + 144(1 - l/y/T+ 5z/3) /625z 2 , with 2 
d 2 ( M ) M //3n(0). 



approximation can be though off as a procedure for eval- 
uating the integral be expanding 1/(1 — gx) in a power 
series in gx and then integrating term by term keeping 
only a subset of contributions. In particular retaining 
the two-point correlations means the following approxi- 
mation, 



dxx 



2n ^— ux 



I 



dxx 



1 



,3/2 ' 



This results in the following result for 1(g), 

i(g)/i(o) 1 



1 - g 2 /uj 



(87) 



(88) 



which has a critical coupling g = g c = y/oJ. The FP 
determinant, however, makes the integral well defined for 
all values of g and thus no critical coupling is expected 
in this case. This is also what happens if Eqs. I)84|) is 
taken into account. As long as fi(0) > the q — value 
of uj(q) does not play a role in determining the position 
of the pole in the FP operator, and we can for simplicity 
assume ui(q) = 0. Then from Eq. I|83|) approximation 
£7(g) = £1(0) for q < (i, we obtain, 



d(q) 



Finally from Eq. 
fi(0) and d(n), 



d(0) 



5 

3 (4tt) 



(89) 



|(HHl we can derive a relation between 



(4tt) ; 



:d 2 (0) 



(90) 



where zq ~ 4 is a root of the nonlinear equation shown 
in Fig. 8 Thus as d(fi) increases so does S1(0) but there 
is no upper limit on d(fi). This is due to the FP de- 
terminant which regularizes functional integrals near the 
Gribov horizon. As J!(0) increases, the transverse-gluon 




10 
k/|i 



FIG. 10: Numerical solution for the vacuum expectation 
value of the FP operator, d(k). The three curves correspond 
do d(k = m) = 38.4 (solid), d(k = fi) = 3.8 (dashed) and 
d(k — fi) — 1.1 (dashed-dotted) respectively, indicating that 
solution may exist for arbitrary choice of d(fi) i.e. no critical 
coupling. 



2-point correlation function decreases at low momentum 
and the ghost corelator function, d(k) increases. This is 
precisely what was found in other gauges using Dyson- 
Schwinger methods and in other approximations to the 
Coulomb gauge. 

In Figs. 9-11 we plot results of numerical solutions to 
the set of coupled equations IM1 1551 I8H 1551 T hese 
should be compared with Figs. 4-6 from Ref. . In 
Fig. 12 we plot the inverse of fi(fc) which is representing 
the transverse-gluon correlation function. As expected it 
is suppressed at low momenta and approaches the per- 
turbative limit as k — > oo. 
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0.60 Lj 




VI. SUMMARY 



10 4 10~ 3 10" 2 10"' 10° 10' 10 z 10' 10* 
k/n 



FIG. 11: Numerical solution for f(k) normalized to f(k 
fi) — 1. Labeling of curves is the same as in Fig. 1101 



10 ' 




In this paper we have studied the role of the Faddeev- 
Popov determinant in the Coulomb gauge. The FP deter- 
minant specifies the measure in the functional integrals 
over gauge field configurations and has so far been ig- 
nored in most calculations of QCD matrix elements in 
the Coulomb gauge. The FP determinant vanishes at 
the boundary of the Gribov region nevertheless it still 
allows for large field configurations near the boundary to 
enhance matrix elements. In particular we have shown 
that the FP operator, corresponding to the running cou- 
pling and the ghost propagator is strongly enhanced in 
the IR, but at the same time no artificial critical coupling 
exists. The same is true for the Coulomb kernel which 
specifies the static, temporal Wilson loop. Finally, the 
instantaneous part of the transverse gluon propagator is 
found to be suppressed as is found in other gauges. 
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10"* L t> <s L, 4, S h h „ 

10 10 10 10 10 10 10 10 10 
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